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Abstract 

The main focus of this research is to study and evalua e the potential of InGaAsSb-AlGaAsSb based 
2pm avalanche photo-detectors. The photo-detector contains a separate absorption and multiplication 
region (SAM) structure. The analysis has mainly been done to understand the electrical response 
characteristics of the devices existing at NASA, and to evaluate alternate structures proposed. 
Calculating the current flow for the existing detector structure, on the basis of its energy band 
diagram, is important. This analysis also helps to find shortcomings in the existing detector structure. 
It is shown that, unfortunately, the existing structure eannot lead to strong multiplication or voltage 
dependent gain. Two alternate structures are suggeste d, that could overcome the inherent flaws, and 
help achieve improved performance. These devices are obtained through modifications of the original 
structure, which include varying the doping levels, and changing the thicknesses of detector sub- 
regions. The results of our study are presented and discussed. 
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1.1 INTRODUCTION 

The main focus of this research is on the potential for 2 |jm avalanche photo-detectors based 
on InGaAsSb-AlGaAsSb, as this material system remains relatively unexplored. The photo-detector 
contains a separate absorption and multiplication region (SAM) structure. The concept is for 
InGaAsSb to act as an absorbing region, and for the AlGaAsSb layer to provide for carrier 
multiplication. The InGaAsSb-AlGaAsSb system is chosen fora variety of reasons. The smaller band 
gap of (E g =0.72) of the InGaAsSb relative to the AlGaAsSb (E g =1.3) region is responsible for 
absorption in that layer, particularly near the 2 pm [1] wavelength. Due to its large band gap, the 
AlGaAsSb region will not absorb incident light at this wavelength. Instead, it should act simply as a 
high-field multiplication region to provide photocurrent gain. The multiplication is based on 
avalanche multiplication upon the application of a reverse bias voltage to the InGaAsSb-AlGaAsSb 
diode structure. The present research is aimed at studying the available structures based on the known 
material parameters and to evaluate device parameters such as dark currents and multiplication, based 
on a physical device model. Model predictions would then be used for comparisons between 
predictions and experimental data. Finally, possible device modifications have also been suggested in 
this report for improving the photodetector response. 

1.2 GalnAsSb PHOTODETECTORS 

As shown recently, it has become possible to use liquid phase heteroepitaxy of lattice matched 
GalnAsSb on GaSb substrates to develop mid-infrared detectors. These quaternary alloy systems can 
provide room-temperature infrared (IR) detectors adjustable for the 1.7-4.5(im range [2]. With 
proper device design, both the cut-on and cut-off wavelengths can be varied over the range of the 
quaternary system. A significant advantage of using a quaternary system over the binary and ternary 
compounds is the ability to vary the band gap while still providing the lattice-matched growth to the 
substrate material. Another expected advantage is that this detector will have good detectivity at room 
temperature, while most commercially available detectors need cooling to liquid nitrogen 
temperatures. 

The GalnAsSb/GaSb structure can be fabricated as either a front side or a backside illuminated 
detector. A diagram of these two designs is shown in figure 1 . The GalnAsSb layers can be grown by 
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liquid phase epitaxy. 




Fig. 1. Front side and backside illuminated detector structures [3], 


1 3 REQUIREMENTS FOR A GOOD PHOTO DETECTOR 

The main requirements of good avalanche photo-detectors are listed below. 

(i) HIGH SIGNAL- 7 O-NOISE RA TIO 

An avalanche photo detector is one which can detects the signals that are coming from a 
source or an object despite their low intensities. The signal coming from the object (for example, 
from a spacecraft) is typically perturbed by the noise. The noise can be introduced during signal 
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propagation, or can also be generated by the detector devices and circuits. In order to detect the 
signal that is largely perturbed by noise, a high signal-to-noise (SNR) ratio detector is desired. A 
high SNR leads to superior detectibility and increases the threshold level for incoming signals. 

(ii) LOW DARK CURRENT 

The dark current of a detector is the current that is flowing when there is no photo-excitation. 
On the other hand, photocurrent is the current when the detector is illuminated by photons. In order 
to detect the photocurrent with precision, the dark current should be small. Typically, the dark current 
should be at least an order of magnitude (or more) smaller than the photocurrent for good detection. 

(in) RESPONSIVITY AT DESIRED WAVELENGTHS 

A good photodetector must have not only a high signal to noise ratio but also superior 
responsivity at desired wavelengths. Responsivity is defined as the ratio of the generated 
photocurrent to the incident signal power. Therefore, the photodetector should yield a maximum 
response for the desired wavelength range and minimum response for all other wavelengths. This 
selective property enables the detector to detect particular signal very effectively and efficiently. In 
the present thesis, the desired wavelength range is around 2 micrometers. 

(iv) LARGE PHOTOCURRENT COMPARED TO DARK CURRENT 

Photocurrent, the current obtained when the detector structure is illuminated, should be larger 
in magnitude in comparison with the dark current for correctly distinguishing the difference between 
them. The photocurrent should be at least an order of magnitude larger than the dark current for 
better detection. 

One way of increasing the detector photocurrent is to engineer internal gain. This can 
typically be achieved by designing the device such that it can undergo carrier multiplication and 
internal avalanche action. Under avalanche conditions, large increases in the current can be achieved 
by very modest increases in the applied voltage. The carrier density typically increases exponentially 
with both the spatial distance and applied electric field. 
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(v) HIGH SPEED OF RESPONSE 

Another useful requirement is the speed of re sponse. Faster speeds imply a larger operating 
frequency, and more data can be processed for digita communication and sensing applications. 

(vi) LOW DOPING LEVELS 

The reduction of doping reduces the dark cu rent as well as the noise in the device because 
of carrier density reductions in the device. This improves the signal to noise ratio and reduces the 
detector dark current. 

(vii) LOW NOISE 

Noise is an unwanted signal that perturbs the detector ability to detect the required signal. 
The amount of thermal noise present is directly propo rtional to temperature. Therefore, noise can be 
reduced by lowering the operating temperature, thereby, suppressing thermal carrier generation. 
Noise is also dependent on the velocity fluctuations, the amount of generation and recombination in 
the system, and doping. Lower doping levels anil reduced defect levels are both helpful in 
suppressing noise. 

1.4 GAlnAsSb AVALANCHE PHOTODIODES 

The structure of aGalnAsSb avalanche photo- detector is shown Fig. 2 [4]. The photo-detector 
consists of n-type InGaAsSb layer that is grown on the n— type Gasb substrate. The InGaAsSb layer 
acts as an absorbing region because of its small band gap. The n-type AlGaAsSb layer that is grown 
on the n-type InGaAsSb layer has a band gap that is large when compared to the InGaAsSb because 
of aluminum composition. This makes it non-absorpti ve particularly at 2.0 pm wavelengths. A p-type 
AlGaAsSb layer is grown on the n-AlGaAsSb layer. The p-n AlGaAsSb layer is responsible for the 
avalanche multiplication upon the application of rev erse bias voltage. All the layers are grown by 
using liquid phase epitaxy and are 3-4 pm thick [5]. 


5 



(^annular contact. 



Fig. 2. The IR SAM-APD structure [4J. 

1.5 ADVANTAGES OF InGaAsSb -AlGaAsSb FOR CURRENT RESEARCH 

The selection of InGaAsSb and AlGaAsSb has several advantages over other materials. The 
potential benefits are stated below: 

1 . The lower band gap of InGaAsSb is responsible for the absorption particularly at 2pm 
wavelengths. 

2. The larger band gap of AlGaAsSb region is responsible for multiplication and avalanching, 
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while serving as an optical window. 

The larger band gap of AlGaAsSb makes it non-absorptive at around 2 pm so the system will 
act as a separate absorption and multiplication region structure. 


2.1 THEORETICAL BACKGROUND 

Before proceeding with the analysis and simulation results of the InGaAsSb-AlGaAsSb 
photodetector system, it is perhaps useful to cover some basic theoretical background. Some of the 
commonly encountered terms and concepts such as the dark current, avalanche multiplication etc., are 
reviewed and presented in the following sections of 1 his report. 

2.2 THEORY OF THE DARK CURRENT 

In this present discussion, a reverse biased p n junction with a depletion width of “W” is 
considered. The electric field “E” can be of arbitrary profile, provided it is sufficiently high in part of 
the depletion layer to cause impact ionization. It is assumed that the electrons and holes traversing the 
depletion layer have ionization probabilities per unit length equal to cl and (3 respectively. These 
ionization probabilities are assumed to be functions of only the electric field E. This means, more 
precisely, that the ionization probability of an electron or a hole is function of its instantaneous 
momentum. A carrier moving in a constant electric field will experience a whole spectrum of 
momentum values, which can be described by a momentum probability distribution for that particular 
value of E. It is only when it has a momentum value towards the high end of the probability 
distribution, that it can be capable of an ionizing collision. The assumption that a or P are a function 
of the local electric field E alone, and not of previou; carrier density relies on the following central 
assumption: The momentum of an individual carrier can sample a large number of values throughout 
the whole momentum probability distribution within a distance that is small compared to: 

1 . The mean distance between the ionizing coui iters, and 

2. The distance in which E changes by an appreciable amount. 

For example, consider an electron, which has just suffered an ionizing collision so that its 
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momentum is reduced to near zero. Before it is capable of creating another electron-hole pair, it 
must gain energy equal to £i, the minimum energy (i.e. the threshold value) required for impact 
ionization. If it gains this energy by moving ballistically through the crystal without collision, 
gaining energy from the electric field, it will require a distance of e, /e*E. Although other 
methods of gaining energy are possible, the distance £j/e*E can be considered as a rough estimate 
of the distance a carrier needs to have its momentum randomized. The requirement #1 (above) 
also implies that the total voltage drop across the region in which appreciable ionization occurs 
must be large compared to £i /e. This condition is assumed to be valid throughout this discussion, 
and generally requires the application of a high device bias. For the sake of simplicity, it is also 
assumed [6] that recombination within the depletion layer can be neglected, although the 
equations to be derived can be generalized easily to include it. 

2.3 DIODE MULTIPLICATION 

Consider a reverse biased avalanche photo-detector in which the high-field region extends 
from x=0 to x=w as shown in figure 3. In this model the electric field E is taken to be such, that 



0 w x 


Fig. 3. Schematic diagram of the high-field region of an APD. 

the electron drift towards the left and holes to the right. Consider, an electron-hole pair generated 
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within the depletion region. 


In traversing distance dx, the electron will suffer on the average ccdx ionizing collisions. 
Similarly, the hole will create an average of pdx hole-electron pairs as it moves a distance dx. These 
secondary pairs will themselves start chains of ionizations, the length of each chain depending on 
where the secondary was created. Thus, if M(x) [7] is the average total number of pairs which are 
generated in the depletion layer as a result of one initial pair being generated at x, then: 

x w 

M(x) = 1+(J a M(x’) dx’) + ( \ PM(x’)dx’), (1) 

0 x 

where oc and p are functions of E which is a function of x\ Differentiating eqn. (1) one gets: 


dM(x) 

= (a - P ) * M(x ) , (2) 

dx 


which has a solution, 


x 

M(x) = M(0) exp[ j (a-P) dx’] . 

0 


w 

= M(w) exp[ J(a-P)dx’]. (3) 

x 


Substituting (3) into eqn. (1 ) for x = w, one obtains 
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w w 

M(w) = 1 - J a exp[- J (a-P) dx’] dx . (4) 

0 x 

Therefore, the multiplication factor M(x) is given by: 

w 

exp[-f (a-p)dx’] 
x 

M(x) = . (5) 

w w 

1 - la exp[-j (a-p)dx’] dx 
0 x 

2.4 TOTAL CURRENT 

In order to calculate the total current, one can consider the rate at which the electron hole pairs 
are generated within an element of distance dx. Since holes and electrons are always generated in 
pairs, only the electrons or the holes (but not both) should be counted. Thus, the pair generation rate 
results in an increase in the hole current at a location x+dx over the current at position x by an 
amount dl p = (dl p /dx) *dx. dl p . This can be expressed as: 

dl p = (al n + pl p + g)dx , (6) 

where al n is the pair generation rate per unit length due to the impact ionization by electrons, pi p the 
rate due to similar collisions by holes, and g is the generation rate per unit length at which the pars 
are generated either thermally or optically. The value for dlp/dx can be obtained from the solution of 
equation (6), and is expressed as: 

x 
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Ip (x)= Ip (0) exp[-J (a-p)dx’] 
0 


x w 

+ J ( al+g ) exp[-J (a-p)dx’] dx , (7) 

0 x 

where I = I p (x)+ I n (x) is the total current and is constant. In quasi-steady-state conditions, the total 
current I is given by: 

W 

1= I p (0)*M(0)+ I n (w)*M(w)+ j g(x) M(x)dx . (8) 

0 

2.5 PHOTO CURRENT 

In the present research, as explained earlier, the AlGaAsSb p-n junction acts as a 
multiplication region and is of particular interest. Assuming that the holes are injected from the left 
side due to photo excitation and no electrons are injected at x=w the above equation is simplified to: 


I = I p (0)*M(0)+ I n (w)*M( w)+ f g(x) M(x)dx , 

0 

W 

= I p (0)*M(0) + 1 g(x) M(x)dx , (9) 

0 

since I n (w) is equal to zero. This assumes that the opi ical signal is incident from the left, and is at an 
energy smaller than the AlGaAsSb bandgap. The current in equation (9) is the total current with 
photo excitation and, therefore, is equal to the photo current. The above expression can be further 
simplified by assuming a equals p. The expression f >r M(x) [8] is then simplified as: 

w 

exp[-J (a-p)dx’] 
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X 

M(x) = 

w w 

1 - la exp [-1 (a-P)dx’] dx 
0 x 


1 

= when ( a=(3 ) . (10) 

w 

1 - | a(x) dx 
0 

The integral in the equation (9) can also be simplified by assuming a constant generation rate g(x). 
The expression for g(x) is given by: 


q*ni*A 

g(x) = , (11) 

Te 

where q is the electron charge, ni is the intrinsic concentration, A is the area of the multiplication 
region and x e is the mean electron life time. Substituting (11) and (10) in (9), the photo-current is 
given by: 

q*ni*A*w 

I = I p (0)*M(0) + . (12) 

w 

V ( 1 - I oc(x) dx ) 

0 

Since M is independent of x when a equals (3, M(0) is same as M(x). Therefore, the above equation 
can be written as: 


I P (0) 


q*ni*A*w 
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I 


(13) 


w 

(1-j a(x) dx ) 
0 


x e *( 1 - 1 oc(x)dx) 
0 


2.6 EXPRESSION FOR THE DARK CURRENT 

When there is no photo excitation, no holes v ill be injected into the detector from the x=0 
side. Therefore, the dark current is the total current when I p (0) is zero, and hence the expression for 
the dark current becomes: 


q*ni*A*w 

Idark = . (14) 

w 

T e *( 1 - 1 0 C(>t)dx) 

0 

The ratio of photocurrent to the dark current is given by: 


I 


photo 


I dark 


i P (0) 

[ 

W 

(1-1 ot(x) dx ) 
0 


q*ni*A*w 

t- ] 

w 

x e *( 1 - 1 ot (x) dx ) 
0 


q*ni*A*w 


(15) 


w 

T e * ( 1 - 1 oa x) dx ) 
0 

By simplifying the above expression, one obtains: 
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Ip(0)*te 


(16) 


I 


photo 


i 


dark 


, for (a = p) . 

q*ni*A*w 


The main aim is usually to increase the above ratio to the maximum extent possible. This can 
be done by decreasing the depletion width “w” and decreasing the intrinsic concentration nj. The 
intrinsic concentration can be reduced by reducing the temperature T because of an inverse 
exponential dependence on T. Reduction of temperature also reduces thermal noise. Though low 
temperature operation can benefit the device performance, it comes at the expense of including low- 
temperature systems. It is thus not a very cost effective procedure, as installing and maintaining the 
cooling equipment requires a lot of energy and capital. 

In general, when a is not equal to p, the ratio of photo current to dark current is given by: 


Iphoto I p (0)*T e *M(0) 


q*A*n;* / M(x)dx 
0 


Since Mavg = l/w*/M(x)dx , 

0 


(17) 


Iphoto 


I dark 


Ip(0)*T e M(0) 

* . 08) 

q*A*n, *w M avg 


Therefore, attempts to maximize the above ratio M(0) requires a condition, in addition to those 
stated above. The electric fields have to be maximum at x = 0 and decrease rapidly with x, since the 
multiplication factor is a monotonically increasing function of the electric field. An electric field that 
rapidly decays from x=0, will ensure that the corresponding multiplication factor will also decrease. 
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Thus, the average value will be much smaller than its magnitude at x — 0. 


2.7 NOISE CONSIDERATIONS 

The noise spectral density within the avak nche photo detector can be calculated by 
considering the above model. Under the assumption that a and P are functions of E only, the current 
generated in the element dx should yield the shot noi ;e, given by: 

< [dip - < dip >] 2 > = 2*e* dl p *df , (19) 

in an element of bandwidth df. This element of noise current can be thought of as originating from a 
noise current generator at x and will undergo the same multiplication as any other current generated 
at x. Thus, if 4> [9] is the noise spectral density ir the total current flowing in the diode, the 
contribution due to the pair generation rate in dx is given by: 

d<Kx) = 2e*M 2 (x) ilp(x) . (20) 


Integrating one obtains: 


v 

<J> = 2e [I p (0)* M 2 (0) + M 2 (w)* I n (w) + q*A*n, * / dlp/dx *M 2 (x) dx] , (21 ) 

0 

where I p (0), I n (w) are the thermally and/or optically generated hole and electron currents 
respectively, entering the two ends of the depletion l ayer. 

The above equation can be further modified by partial integration of the last term. This 

yields: 


w w 

i dlp/dx *M 2 (x)dx = - I p (0)* M 2 (0)+ M 2 (w)* I p (w) - 2 / dlp/dx *M 2 (x)dx. (22) 
0 0 
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By rearranging the above equation one can get the final expression for the spectral density. It is 
given by: 

W <J) 

= I* 2*e*[ 2 */ a M 2 (x) dx - M 2 (w)] + 4*e* I p (0)* M 2 (0) + M 2 (w)* I n (w) + 

0 

w 

+ 4’*‘e 1 " | g *M 2 (x)dx , (23) 

0 

where q is the electron charge, l p (0)and I n (w) are the hole and electron currents entering the 
depletion region at x=0 and x=w respectively, g(x) is the rate at which the hole-electron pairs are 
generated (thermally and/or optically) in the depletion region, and M(x) is the average gain associated 
with a hole electron pair generated at position x. Provided the values of cc(E), (3(E), and E(x) are 
known, the above equation can be used to compute the noise spectral density for any desired 
distribution of injected carriers. 

CASE #1: 

We consider the case when a is equal to (3. When there is no photo excitation l p (0)and I n (w) are zero. 
Assuming the g(x) factor [10] is due to thermal excitation, the above equation reduces to: 

w w w 

4>da = 2*e*[ 2*1 a M 2 (x)dx) ] + j g *M 2 (x)dx *[2* ja M 2 (x)dx- M 2 (w)] , (24) 

0 0 0 

where <t>da is the noise spectral density associated with the dark current generated in the region 0 < 
x < w . The primary dark current will then be given by: 

w 

Ida = 1 g (x)dx . (25) 

0 
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The mean dark current at the output is given by: 


w 

(Ida)out = J g(x)M(x)dx . (26) 

0 

The distributed dark current can be replaced by an equivalent photocurrent source in the p + 
region with the same mean value Ida . This can be done by introducing an effective average gain M e fr 
, and an effective noise factor F e ff , differing from those for the injected currents, chosen such that 
Ida lout and <|)da are unchanged. Hence, M e ff would be given as: 


(Ida)oi t 

M e ff — 

Ida 


(27) 


The above equation can also be written as : 


\ g(x)M(x)dx 
0 

Me ff = . (28) 

\* 

! g (x)dx . 

0 


Similarly, the F e ff factor would be given by : 


F e fr 


i|>da 


2 *Cj* Ida* M 2 eff 


(29) 
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w 


w 


w 


w 


[ fg (x)dx ] [2* Jg (x)M 2 (x)dx +[ fg (x)M(x)dx ] *[ fa* M 2 (x) *g (x)dx -M 2 (w) ]] 
0 0 0 0 


Fef f - (30) 

w 

[ f g (x) *M(x)dx ] 2 
0 

By using the relation a =P, the F e ff factor turns out to be as: 

F eff = M eff (31) 

Consequently, the noise spectral density would then be given by: 

<t> = 2*q*I da *M 2 eff • (32) 


CASE #2: 


If the electron and hole ionization coefficients are not equal, but vary in roughly the same 
manner with the electric field, then the approximation p = ka is reasonable one, with k some fixed 
constant. With this approximation, the integral in the equation (24) can be performed [11]: 

w w 


2 * f a M 2 (x)dx = 2/( 1 -k) *J (a-p) M 2 (x) dx 
0 0 

M 2 (w)- M 2 (0) 


1-k 


(33) 


Therefore, equation (24) becomes: 


<(> = I* 2*e*[k M 2 (w)- M 2 (0)] + 4*e* I p (0)* M 2 (0) + M 2 (w)* I n (w) + 


18 



1-k 


w 

+ 4*e*Jg* M 2 (x) dx . (34) 

0 

In the usual case of interest, most of the ionizaion occurs over a very short distance that may 
be only a small fraction of the total depletion width. If a and b represent the locations in the depletion 
layer of the two ends of the region in which the most ionization occurs, one can write M(x) =M(a) for 
0 < x < a, and M(x) * M(b) for b < x < w. Therefore, the above equation can be written as: 

I* 2*e*[k M 2 (b)- M 2 (a)] + 4*e* I p (a)* M 2 (a) + M 2 (b)* I n (b) 

♦ = 

1-k 

w 

+ 4*e* J g *M 2 (x)dx , (35) 

0 

w 

where I p (a)= I p (0) + J g(x)dx , and 

0 

w 

In(b)= ln(w) + J g(x )dx . (36) 

b 

The above equation for noise spectral density can be simplified as: 

<j> = M 3 (a)*I p (a)* [1 + ( 1 -k) *( M(a)- l) 2 /k*M 2 (a)] + 

+ M 3 (a)*I p (a)*[l -(l-k)*( M(b)- 1) 2 /M 2 (b)] 
b 

+ 4*e* J g *M 2 (x)dx 
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a 


w 


+ k M 2 (b) - M 2 (a) * l g *M (x)dx . 
0 


1-k 


( 37 ) 


Here the relationship that M(a)-1 = k( M(b) -1 ) is used , which is true in the case of p = ka. 
CASE #3: 

We now consider the situation where P = ka and M °=. For large values of M(b) and M(a) , 

the M(a) =k*M(b) and the last term in the equation is the dominant one. Therefore, <j> is given by: 

<J> = 2*e*I*M(a)*M(b) , M(a),M(b) ►«. (38) 

This is the spectral density of the current noise generated in the depletion region. 

In summary, in the case of a = p, the noise spectral density is proportional to the third power 
of the diode multiplication. Since the photodiode signal will be proportional to M, the square of the 

O'! 

signal-to-noise ratio will be proportional to M /(2*e*I t *M +<J>o ) [12] where <|)o is the equivalent 
spectral density, referred to the input, of all other noise sources. This will have a maximum for M 
=2*<|>o /2*e*I t . If ka is equal to P, roughly the same comments apply. However, the noise spectral 
density can be greater or less than the 2*e*I t *M 3 ’ depending on where the light is absorbed in the 
junction. If it is absorbed mostly to the left of the high-field region, then the injected current I ln « I t 
and: 


4> / 2*e*Ijn =M 3 [l+l-k/k*((M-l)/M) 2 ] . (39) 

Conversely, if most of I t originates from light absorbed to the right of the high field region, 
then I ^ * I t and: 


4> / 2*e*I in =M 3 [l-(l-k)*((M-l)/M) 2 ] . 


(40) 
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Equations are plotted in Figs. 4-5 for various values of k. From the graphs, the following 
conclusions can be drawn: 



Fig. 4. Noise spectral density/2eli n versus Multiplication for either injected holes 
or electrons if P = Ka . 
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Fig. 5. Noise spectral density/2eli n M 3 versus multiplication for either 
injected holes or electrons if p = Ka . 


1 . The noise at a given value of M will be less if it can be arranged that most of the carriers 
entering the high field region are of the more strongly ionizing type. 

2. If one wishes to compare the noise measurements made on a real diode with the theoretical 
curves for (3 = ka, k should be considered an adjustable parameter roughly equal to the (3(E 
max)/ aCEmax), where En, ax is the maximum field in the depletion layer for a given value of M. 
However, the figure shows that for the more strongly ionizing carrier, F is fairly insensitive to 
the value of k for low to intermediate values of M. Hence, the theoretical curve for F vs k for 
large M may be indistinguishable from a plot of experimentally measured values of F vs M. 
For the less strongly ionizing carrier, on the other hand, F is fairly sensitive to the value of k, 
so that a single value of k is not expected to give good agreement with experiment, except 
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when the approximation (3-ka happens to hold. 


Measurements of the hole impact ionization rate in Al x Gai. x Sb for different compositions [13] 
is shown in the figure. The fractional alloy composition of As in the Al x Gai. x AsySbi. y is 0.08 (i.e. 
y = 0.08) for the photo-detector structures of interest to NASA. Since the relative composition of 
As is low, it can be neglected, and therefore, the properties of Al x Gai. x Sb can nearly be attributed 
to the actual AlxGa^AsySb,^ . For the given detector, the composition (x) is close to 0.26, for 
which the k=(f}/a) ratio is approximately equal to 2.1 . 


Alloy Composition (x) 
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Fig. 6. The variation of k=(|3/a) with the composition of Al x Gai. x Sb [13]. 


3.1 ELECTRIC FIELD CALCULATIONS 

Neglecting the GaSb substrate of the photo-dotector, the device contains only three regions 
namely, the n-type InGaAsSb layer, the n-type AlGAAsSb and p-type AlGaAsSb layers as shown in 
Fig. 7. The energy band diagram (Fig. 8) and the elect ric field at each point can be calculated at each 
point by assuming a certain electric field and using Poisson’s equation at each junction. 
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Fig. 8. The expected energy band diagram of the above detector. 
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The electric field at x=0 + is assumed to be a value equal to EO. The electric field can be positive or 
negative. Here, two cases have been considered, one with a positive electric field and another with 
negative electric field. The equation governing the tv o junctions is expressed as: 

-8el-8v 1+0.75 -8v2-Eg2- 8e2 =Va, (41) 

where 8el is the energy gap between the conduction band and Fermi level in InGaAsSb and 6e2 is 
the energy difference between the valence band and Fermi level in p-AlGaAsSb. The parameters 8v 1 
and 8v2 are the voltages developed in the two junctions, while 0.75 represents the band offset 
between n-InGaAsSb and n-AlGaAsSb, and Eg 2 is ‘he band gap of AlGaAsSb. This equation is 
derived from the continuity of energy across the entire device, and is valid for a given applied voltage 
Va. 


CASE #1: The procedure for evaluation of the electric field profile and voltage 
distribution within the device is described in this section. This procedure has been used throughout 
this thesis for calculations of the energy band diagrams. For concreteness, we begin by assuming that 
the initial electric field EO is greater than zero, and the voltage 8vl is negative. Use is made of the 
governing equation is the Poisson equation. Not only EO, but also a negative value of 8vl, is 
assumed. By applying Poisson’s equation for x less lhan zero: 

dE/dx q*[Nd InGaAsSb ' " n(x)]/ £ InGaAsSb > (42) 

where n(x) is the carrier density for x less than zero. The electric field is varied according to the 
above equation with respect to x by incrementing x in small steps with a starting value at x=0. The 
initial value of the electric field can be obtained from the continuity of D (the electric flux density): 

E(O’) = £ InGaAsSb* EO £ AlGaAsSb . (43) 

The initial value of carrier density can be calculated from the formula: 
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n(0)=ndl *exp((q*8vl)/(kb*T)) . 


(44) 


The electric field is then incremented by using the Poisson’s equation: 

E(x-8x) = E(x)- q*[Nd foGaAsSb - n(x)]* [8x]/ e MiaAsSb • (45) 

The voltage in the elementary length 8x is given by: 

V(-5x) = [E(x)+ E(x-8x)]* [8x] 12 . (46) 


In order to proceed, the carrier density n(x) is incremented by using the following equation: 

n(x-8x )=n(x)*exp((-q*V(-8x)/(kb*T)) . (47) 

Poisson’s equation is solved until the electric field becomes zero at x = -dl , which is the distance of 
the depletion region of interest. The differential voltages obtained are added to get a final voltage V. 
This voltage is compared with the assumed value of 5vl. If both the voltages have the same 
magnitude and opposite sign then the chosen EO is correct. The difference in the sign of 8vl and V is 
attributed to that the V is measured in reverse order compared to 8vl . The chosen EO is varied until 
both the voltages match. 

Once the 8vl is matched by adjusting EO, then the calculations are carried out in the 
AlGaAsSb region, and Poisson’s equation is used once again as : 


dE/dx = q * [Nd AlGaAsSb - n(x)]/ e AlGaAsSb- (48) 


The adjusted EO is now used as the initial value for this equation. The initial value of carrier 
density is given by the following equation: 
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n(0) = ncl*exp((q*(8vl +8» l-0.75))/(kb*T)) , 


(49) 


where ncl is the carrier density in the conduction band of InGaAsSb. Poisson equation is solved in 
the same way as explained above. The electric field if. incremented by using. 

E(x+8x)=E(x)+([8x]*q*(nd 2-n(x))/( £ AiGaAsSb) • (50) 

The voltage in the elementary length 8x is given by: 

V(8x) = -[E(x)+E<x+8x)]*[8x]/2. (51) 

The carrier density is then modified according to: 


n(x+5x )=n(x)*exp((q*V(8x)/(kb*T)). (52) 


Poisson’s equation is solved until the distance x becomes x = 1, where 1 is the given thickness of 
the AiGaAsSb region. The voltage V obtained is compared with the 8v2 obtained from the 
equation: 


-8el-8vl+0.75 -SvS-Eg- &2 =Va. (53) 


If V matches 8v2, then the chosen 8vl is deemed to b< : correct and E0 is the accurate electric field. If 
not, a different 8vl is chosen and the whole procedu re is repeated until the all the values match. 

CASE #2: When E0 is less than zero , the same above procedure is followed except 

that the 8vl is chosen to be positive. The above procedure can be explained by the following 
flowchart of Fig. 9. 
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Flow chart for calculating the electric field, 








5.1 RESULTS AND DICUSSION 

A.summary of the various results obtained in s he analysis of the detector is provided in this 
section. The analysis has mainly been done to unden tand the electrical response characteristics of 
the existing structure, and to evaluate alternate structures. Calculating the current flow for the 
existing detector structure, on the basis of its energy b, ind diagram, carries out the response analysis. 
This analysis also helps to find imperfections and shoitcomings in the existing detector. Once such 
drawbacks are recognized, these can be overcome through the design of alternate geometries. It is 
shown that, unfortunately, the existing structure at NA SA is lacking in these properties. Two alternate 
structures have been suggested, that should overcome the flaws inherent in the original detector 
design, and help achieve improved performance. These structures are obtained through various 
modifications of the original structure, which include varying the doping levels, and changing the 
thicknesses of detector sub-regions. 


5.2 ENERGY BAND DIAGRAM OF THE EXISTING NASA DETECTORS 

The energy band diagram for the existing detector can be obtained by combining the energy 
band diagrams of the three individual junctions of the detector. The junctions are n-GaSb:n- 
InGaAsSb, n-InGaAsSb:n-AlGaAsSb, and n-AlGaA ;Sb:p-AlGaAsSb. For convenience, these are 
referred to as the first, second and third junctions, respectively. 
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Fig. 10. The Energy band diagram of GaSb-InGaAsSb junction. 

Figure 1 0 shows the partial energy band diagram (i.e. only the conduction band) for the n- 
GaSb:n-InGaAsSb junction. The valance band has been omitted for simplicity, since it lies E g below 
the conduction band. The shape of the energy band diagram is typical of n-n heterojunctions [14, 
15]. The energy variation is maximum at the origin (distance =0) and is equal to the band offset, 
which has a 0.325 eV value. 
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Fig. 1 1. The energy band diagram of InGaAsSb-AlGaAsSb junction. 


The energy band diagram (conduction band' for the second junction is drawn in the same 
way as the first junction. Note that the origin of the second junction is at a distance d=3 (am from the 
origin of first junction (i.e. d = 0). This corresponds to the thickness [1 6] of the n-InGaAsSb region, 
which, for this detector is 3 pm. 

The energy band diagram for the third junction resembles that of the conventional p-n 
junction. The origin of this diagram is at a distance d r 6 pm, since the AlGaAsSb [17] region in the 
NASA device is 3pm thick. Here, the energy band diagram penetrates more into the n-AlGaAsSb 
region as compared to the p-AlGaAsSb, because the p-AlGaAsSb region has a larger doping with 
respect to compared to n-AlGaAsSb region. Note that all the three energy band diagrams have been 
drawn with no external voltage applied i.e. when applied voltage is zero. 





Fig. 12. The energy band diagram for AlGaAsSb p-n junction. 

Once the energy band diagram for each junction has been obtained, the full energy band 
diagram can be deduced by combing each of the three individual segments. The combination has to 
be based on the condition that, when the current due to generation rate Go is neglected, the current 
flowing in each junction should be same for continuity. Consequently, the applied voltage will then 
divide between the junctions in such a way that the current at any point is the same and is equal to the 
total current flowing in the device. The idea similar to putting together a collection of resistors 
connected in series, subject to a constant total current, and a voltage division between the constituent 
resistive elements. 
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Fig. 13. The combined energy band diagram for the first two junctions. 


Figure 1 3 shows the energy band diagram ft >r the first two junctions based on the above 
voltage dividing principle. The voltage across the twc junctions is divided in such a way that the total 
current is constant and is equal to 1 .358 *10 3 A/m 2 . 1 his voltage division can be done from the J-V 
characteristics of the two junctions. From this, one obtains: 

Total applied voltage: V = V] + V 2 - 0.30009171 V , (54) 


and, total current: J « Ji = J 2 ® 1 .358 *10 3 A/m 2 . (55) 












Fig. 14. The complete energy band diagram for the detector at V = 0 V. 
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Fig. 15. The complete energy band diagram for the detector at V - 0.4 V. 

Combining the energy band diagrams, as shown figure 14, draws the complete energy band 
diagram of the detector when no external voltage applied. The energy band diagram for the 
detector when an external, non-zero bias is applied is shown in figure 15. The applied voltage is 
equal to 0.400641 V. Most of the voltage is predicted to drop only across the third junction (V3 = 
0.4V), to maintain a constant current equal to 5.1) x 10 -4 A/m 2 in the detector. The negligible 
voltage drop in the first and second junctions can be explained on the basis of the J-V 
characteristics of the first and second junctions. 

5.3 COMPAR1SION OF SIMULATION Ol JTPUT WITH NASA DATA 

In this section, results of photodt tector responsitivity for the existing NASA 
devices are compared with theoretical curves obtained from our simulations. The responsivity 
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measurements for the original NASA structure were obtained at Astro Power Inc. and other 
laboratories. Illumination on the detector could be incident from either the front side or the back 
surface. By convention, for front side illumination, light entering the detector would fall on the p- 
AlGaAsSb. Conversely, backside illuminated would be the situation for the light beam to be 
incident from the GaSb-side. The responsivity of the detector obtained from experimental 
measurements with front- and back-side illuminations, is shown in Fig. 16. 



Fig. 16. Spectral response of front side and backside illuminated NASA structure 

[18]. 

As seen from the figure, the back-illumination curve has a maximum responsivity 
around the 2000 nanometers (nm) wavelength. The responsivity falls off quickly as the wavelength 
changes from 2000 nm. The “front” curve, on the other hand, has a responsivity that is constant 
from 2000- 1 700 nm, and then the responsivity begin to decrease slowly for lower wavelengths. The 
difference in their responsivities is due to the material that is primarily responsible of the photon 
absorption. For “backside” illumination, GaSb absorbs most of the light, particularly below the 
2000 nm regime. This causes the resonsivity to fall sharply below 2000 nm in this mode of 
operation. 
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The absorption coefficient (a) of InGaAsSb for these curves is given to be 5.5 *10 5 m' 1 at 
2000 nm. However, a is not a constant and chang es with wavelength due to absorption by the 
band-edge “Urbach” tails, and increases in the density of states with energy. The variations of a 
were not measured, and neither are there any reports in the published literature. Given this 
unknown, the wavelength dependent alpha values were calculated from the responsivity of the 
“frontal-illumination” data curve. The alpha was asi-umed to be nearly constant in the region below 
2000 nm. The resulting alpha variation thus obtained from the responsivity measurements, is shown 
in Fig. 17. 



Fig. 17. Alpha versus wavelength for the existing NASA detector. 
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Fig. 18. Comparison of experimental responsivity with the existing data. 

Shown in Fig. 18 are the experimental data for the normalized responsivity, and a 
theoretical calculation based on the calculated alpha curve of Fig. 1 7. As can be seen, the 
theoretical predictions match both the front and back curves fairly well. Deviations at the long 
wavelengths are probably due to device leakage currents. 


38 






Fig. 19. Comparison of experimental and NASA curves for front illumination. 


For alpha assumed to be nearly constant foi wavelengths below 2000 nra, the theoretical 
responsivity obtained is shown in Fig. 19. The experimental data roughly matches the calculated 
curve. The experimental curve is also drawn for the backside illumination case, for completeness. 
As seen from the figure, the experimental curves ar; nearly identical for both the front- and back- 
side illumination situations for wavelengths greater than 2000 nm, while they differ for lower 
wavelengths. This difference is arises because upon backside illumination of the device, the 
absorption from the GaSb layer causes an additional decrease in photon intensity. This causes the 
responsivity to fall more quickly. The experime ntal back curve has the same shape as the 
theoretical curve, but falls more rapidly after 1800 nm wavelengths. 
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The other structures available and tested at NASA were also studied, and their responsivity 
values compared. For example, wavelength-dependent responsivity for NASA structure D19-K204- 
1 1#8 , had a response as shown in Fig. 21 . 
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Fig. 21. Responsivity versus wavelength for D19-K204-ll#8 detector. 


The responsivity versus wavelength curves for the NASA detector D19-K204-1 1#8 at the 
two applied voltages of OV and IV, are shown in figure 21. As seen from the plots, an increase in 
voltage causes a corresponding increase in the responsivity for a given wavelength. These 
experimental data are compared with the theoretical predictions obtained from the extracted alpha 
values of figure 17. The resulting plots are shovti in figure 22. As seen from figure 22, the 
experimental data and calculated responsivity match almost exactly for V = OV, and also are in 
nearly agreement at the IV bias. This validates the alpha values and bears out the theoretical 
modeling. 
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Fig. 22. Comparison between experimental and existing curves for D19-K204-1 1#8. 


5.4 ANALYSIS OF EXISTING DETECTOR STRUCTURE 

The characteristics, including the energy band diagram, for the photo-detector existing at 
NASA LaRC have been explained in the previous section. The detector was designed for several 
NASA applications mainly at the 2pm wavelength. In order to enhance the sensitivity, it is highly 
desirable to achieve internal carrier multiplication and avalanche action for amplifying the photo- 
current. Unfortunately, the current detector described in the previous section seems to have several 
flaws in its design. First, the dark current of the detector, which should ideally be small for low 
background noise detection, is actually large. The dark current in this detector is in the same orders of 
magnitude as the photo current. This makes it very difficult to distinguish between the actual 
photocurrent response, and the dark currents. The detector also suffers from the problem of being 
able to achieve high carrier multiplication. The avalanching property is useful because, in most of the 
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detectors, the incoming signal will be weak. This is especially so if the incoming photon signal is 
weak due to atmospheric attenuation. For good detect on, such as weak signal needs to be amplified 
to produce a strong, distinguishable output. Therefore, without avalanche action and internal gain 
through carrier multiplication, it is very difficult to delect weak signals. This is, therefore, is one very 
important drawback of this detector that needs to be corrected. The following section present 
alternate structures that attempt to rectify the above faws. 

5.4.1 GUIDING PRINCIPLES FOR GOO D PHOTO-DETECTION 

The following steps should be taken for improving the performance of the photo-detector. 
The main aim is to reduce dark current and increase t he photo-current. 

1. Reduce the doping N d within the n-InGaAsSb region. The reduction of doping N d will 
increase the width of depletion layer in that region. This would also cause a decrease in the 
carrier recombination, since relatively fewer free carriers would be available for direct 
recombination. The reduction in carriers als< > brings about a reduction in the dark current. 
This would effectively work to suppress the background noise. This is the main advantage of 
reducing the doping in InGaAsSb region. Reducing doping also increases the resistance in 
that region which will increase the IR voltage drop. Consequently, the electric field would 
increase. This in turn, can facilitate an increase of the drift current, and the possibility for 
carrier impact ionization. 

2. Reductions in doping within the n-AlGaAs >b region, and to make it nearly intrinsic, if 
possible. This reduction in doping would help; in reducing 6Ev [22] between the InGaAsSb- 
AlGaAsSb regions by pushing up the AlGaAsSb energy band diagram. This will also help 
increasing the electric field necessary to drive the holes and to cause impact ionization in the 
p-n junction. Finally, the possibility of tunneling in p-n junction can also be avoided by 
reducing the doping. 

3. Reducing the width of AlGaAsSb region. This reduction will cause the electric field in the n- 
p junction to penetrate deeper into the InGaAsSb region. This also reduces possible 
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recombination. Finally, the reduced width will lower the carrier transit time, thereby, 
increasing the frequency response and speed. 


5.5 ANALYSIS OF TWO ALTERNATE DEVICE STRUCTURES 

In keeping with the above strategies, two alternate photo-detector structures are proposed for 
improving the performance of the existing device. In both the alternate structures, the n-GaSb 
region is completely removed, as it really does not serve any useful role in either the carrier 
absorption, or multiplication processes. The two alternate structures are: 

1 . An intrinsic InGaAsSb-Intrinsic AlGaAsSb- p AlGaAsSb detector, and 

2. An intrinsic InGaAsSb-p (lightly doped and almost intrinsic) AlGaAsSb-p AlGaAsSb 
detector. 

The analysis of the above two structures is done, in an attempt, to see if the performance 
can be improved. The first structure consists of very lightly n-InGaAsSb (almost intrinsic) with a 
doping concentration of 1 0 22 m' 2 , and n- AlGaAsSb (almost intrinsic) with a doping concentration 
of 5 x 10 21 m' 2 , and a p-type AlGaAsSb region. Note that the width of the AlGaAsSb region is 
reduced to 1 pm from the 3 pm dimension of the original detector. 
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Fig. 24. The electric field for the modified structurel at V - OV. 
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The energy band diagrams of the above structures were calculated at different voltages. The 
energy band diagram with no applied voltage is shown in figure 23. Note that the energy band 
diagrams of n-InGaAsSb-n-AlGaAsSb and AlGaAsSb p-n junction are combined, and there is a 
definite inter-penetration of electric field from the AlGaAsSb p-n junction into the InGaAsSb region. 
The corresponding internal electric field profile is shown in Fig. 24. The discontinuity at the 
InGaAsSb-AsGaAsSb junction is due to the difference in dielectric constants between the two 
materials. The energy band diagrams for other applied voltages ranging from 0.5 V to 7V have also 
been shown in Figs. 25-29. As seen from the figures, the increase in the voltage causes the energy 
band diagram to shift up in InGaAsSb-AlGaAsSb region. For voltages below 2 V, there is barrier to 
hole flow in the n-AslGaAsSb region. Hence, the photocurrents at such biasing levels are not 
expected to be very large. However, as the biasing voltage reaches beyond a 2V value, the valance 
band barrier disappears. Hence, it should be possible for holes in the valance band of InGaAsSb to 
move unimpeded, and be transported to the right, thereby contributing to a net current. The very high 
electric field will promote high drift velocities, and could even cause impact ionization at the 
AlGaAsSb p-n junction. The width reduction of AlGaAsSb will enhance this, not only because of the 
higher electric fields, but also of negligible recombination over a narrower width. 
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Fig. 28. The complete energy band diagram of the modified structure! at V = 2V. 
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Fig. 29. Complete energy band diagram of the modified structurel at V = 7V. 


The above structure is also predicted to havt a low dark current because of reducing doping. 
This bodes well for the overall selectivity and detect ivity of the device. Also, a large photocurrent 
should result, given the possibility for carrier avalanche at high biasing, and low recombination along 
the transit path. The large band bending is indicative of the potential for high impact ionization. 





Fig. 30. The complete energy band diagram of the modified structure2 at V = 0V. 

A second structure is now proposed and investigated. The second structure consists of very 
lightly n-InGaAsSb (almost intrinsic) with a doping concentration of 1 0 22 m' 3 , and p-AlGaAsSb 
(almost intrinsic) with a doping concentration of 5 x 1 0 21 m' 3 , and a p-type AlGaAsSb region. Note 
that, similar to the first structure, the width of the AlGaAsSb region is maintained at 1pm. The energy 
band diagrams for the second structure were also calculated, and the shown for various applied 
voltages in the following diagrams. Fig. 30 is for a zero bias case, the diagram of Fig. 31 for 0.5 V 
and Fig. 32 for 2 V bias. 
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In this second structure, as the applied voltage increases, most of the voltage is predicted to 
drop across the AlGaAsSb p-n junction. Like the first structure, this device scheme structure would 
also help drift holes in the InGaAsSb and AlGaAsSb valance bands. There is also a distinct 
possibility to cause impact ionization across the AlGaAsSb junction, and to produce a large 
photocurrent. 

Since the success of these photo-detectors depends on the degree and extent of impact 
ionization for internal gain, it is necessary and germane to evaluate the ionization rates. The hole and 
electron impact ionization rates for the AlGaAsSb region of the existing NASA detector were, 
therefore, calculated for different applied electric fields and for different compositions. As discussied 
previously and shown in Fig. 6, the change in the composition, has no effect on the hole impact 
ionization coefficient (J3) [23-25]. However, the electron 



Fig. 33. The hole impact ionization rate versus inverse of electric field. 
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Fig. 36. The electron impact ionization rate versus inverse of electric field. 

ionization coefficient (a) changes by a factor k (p/a), and equals 2.1 for the given 
composition of A1 in AlGaAsSb for the existing NASA detector. 

Figs. 33 and 34 show the calculations for electric field vs ionization coefficient for 
holes. The electric field dependence of the electron ionization coefficient a is shown in Figs. 
35-36. Both the E vs a and 1/E vs a graphs [26] are plotted for convenience. Since the 
impact coefficient varies with alloy composition, a set of possible A1 compositions in 
AlGaAsSb, has been used. As seen from the graphs, the increase in composition increases the 
electron impact ionization rate, for the same applied electric field. 

Calculations of multiplication factor [27] for the AlGaAsSb diode are shown in the 
following graphs. As seen from Figs. 37-40 the multiplication factor, M, increases with 
increasing A1 composition x in the AlGaAsSb material for a given applied voltage. The 
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multiplication factor also increases with the applied voltage for the fixed composition. These 
results clearly show that the alloy composition is also an important parameter, and should be 
carefully considered in the design of a good avalanche photo-detector for this 2 pm 
application. 






, - • - , r-, , . „ . 

Fig. 39. The variation of multiplication factor with alloy composition at V = IV. 
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